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Abstract. We show that a helical shear flow of a magnetized plasma may serve as an efficient amplifier of Alfven 
waves. We find that even when the flow is purely ejectional (i.e., when no rotation is present) Alfven waves are 
amplifled through the transient, shear-induced, algebraic amplification process. Series of transient amplifications, 
taking place sequentially along the flow, may result in a cascade amplification of these waves. However, when a 
flow is swirling or helical (i.e., some rotation is imposed on the plasma motion), Alfven waves become subject 
to new, much more powerful shear instabilities. In this case, depending on the type of differential rotation, both 
usual and parametric instabilities may appear. We claim that these phenomena may lead to the generation of 
large amplitude Alfven waves and the mechanism may account for the appearance of such waves in the solar 
atmosphere, in accretion-ejection flows and in accretion columns. These processes may also serve as an important 
initial (linear and nonmodal) phase in the ultimate subcritical transition to MHD Alfvenic turbulence in various 
kinds of astrophysical shear flows. 
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1. Introduction 

Most of astrophysical objects involve different kinds of 
plasma flows. Recently it was fully realized that collec- 
tive phenomena in flows with spatially inhomogeneous ve- 
locities (shear flows hereafter referred as SF) are charac- 
terized by remarkable, so called "nonmodal" processes, 
related with non-self-adjointness of linear dynamics of 
perturbations. Namely, it was found that SF: exchange 
energy with sound waves (Butler & Farrell 1992); couple 
different collective modes with one another and lead to 
their mutual transformations ( phagelishvili ct al. 1996D ; 
generate nonperiodic, vortical modes of collective be- 



mas 



dRogava et al. 1998 

Mikhailenko et al. 2000D " 



Volponi ct al. 2000 



magnetized plasmas 



haviour (so called "shear vortices" Rogava et al. 1998) 
- which eventually may or may not acquire wave- 
like features; excite beat wave phenomena both in 



neutral fluids ( Rogava fc Mahajan 1997 ) and plasmas 



(Poedts et al. 199S 



These processes take place not only in neutral fluids 



(Butler & Farrell 1992 




Rogava & Mahajan 1997 


), Stan- 


dard MHD (|ChagcUshviU et al. 1996, 


Rogava et al. 1996, 


[Tatsuno et al. 2001) and electrostatic plas- 
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but also in strongly 
with anisotropic pressure 
( Chagelishvili ct al. 1997), electron-positron plas- 
mas (Mahajan et al. 1997) and dusty plasmas 
(Poedts et al. 200C). The possible role of these phe- 
nomena in astrophysical context was immediately 
realized and a number of astrophysical applica- 
t ions, including pulsar magnetospheric plasmas 
( Mahajan et al. 1997 ) , solar atmosphe ric phenomena 
( Poedts et al. 1998 , Rogava et al. 2000 ) an d galactic 
gaseous disk dynamics ( Rogava et al. 1999| ) appeared 
within a span of years. 

One of the main shortcomings of all these studies, 
stemming from the very design of nonmodal schemes, is 
that the description is given in the wave number space 
(k— space), the knowledge about the appearance of shear- 
induced phenomena in the real, physical space is lacking. 
Two more serious limitations of these investigations are 
related with the neglect of the back reaction of perturba- 
tions on the mean flow and with the omission of viscous 
dissipation effects. However recently a successful effort was 
made in the direction of spatial visualization of the shear- 



induced wave transformations (Bodo et al. 2001), became 
clear that shear-induced processes in SF are quite robust 
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and easily recognizable, even in the presence of quite heavy 
dissipation. 

All these processes were studied predominantly for 
simple, plane-parallel flow geometries with linear veloc- 
ity shear profiles. However, recently, a new method was 
developed ( Mahajan fc Rogava 1999| ), which allows a lo- 
cal analysis of the dynamics of linearized perturbations 
in SF with arbitrarily complex geometry and kinematics 
'U{x,y,z). It was found that a slight deviation from the 
plane-parallel mode of motion brings into the game a va- 
riety of new, exotic and asymptotically persistent modes 
of collective behavior. Still the systematic investigation of 
shear-induced processes in kinematically complex SF is a 
challenging task in a state of infancy. The original study 
has been limited to two-dimensional flow patterns of neu- 
tral fluids ( Mahajan fc Rogava 1999 ). 

One begins to wonder whether those astrophysical 
systems, where kinematically complex modes of plasma 
motion are definitely present — e.g., astrophysical jets 
( Ferrari 1998|) or s olar tornados ( Pike fc Mason 1998 , 
Velli fc Licwer 1999 ) — might sustain these collective phe- 



nomena and what kind of observational consequences 
could they lead to. 

In this paper we consider interaction of a helical flow 
with Alfven waves generated within the flow. This prob- 
lem is actual in plasma physics and plasma astrophysics 
from experimental/observational, theoretical and numeri- 
cal points of view. 

We are used to observe Alfven waves, being in inter- 
action with plasma flows, in different astrophysical situa- 
tions. For example, for the long time it is known that the 
Sun radiates Alfven waves: outward propagating Alfven 
waves are routinely observed in the solar wind flow at 
r > 0.3 AU ( |Hollwcg 199C| ). The observations of quasiperi- 
odic pulsed ionospheric flows (PIF) have shown that the 
PIF are driven and correlated with Alfvenic fluctuation s 
observed in the upstream solar wind ( Prikryl et al. 200^ . 
While the observations of the solar transition region 
( Peter 2001 ) imply that structuring of the transition re- 
gion involves closed loops and coronal funnels showing 
unambiguous evidence for the presence of passing Alfven 
waves. 

In the theoretical domain the problem of the inter- 
action between Alfven waves and ambient flows is quite 
popular topic of studies in the wide range of applications 
including different kinds of laboratory, geophysical and as- 
trophysical plasma flows. Recent studies of toroidal flows 
in axisymmetric tokamaks ( van der Hoist et al. 200C ), for 
instance, revealed that these flows generate low-frequency 
Alfven waves. The problem of Alfven waves sustained 
by plasma flows is especially popular in the context 
of solar physics. It is argued that chromosphere and 
transition region flows are primary energy sources for 
the fast solar wind; while Alfven waves, generated by 
post-reconnection processes, are continuously interacting 



/3 dBalbus fc Hawley 1991D and low-/? ( [Tagger et al. 1992| ) 
disks. Most recently it was claimed that the Accretion- 
Ejection Instabihty ( [Tagger fc Pellat 1999| ) can extract 
accretion energy and angular momentum from magnetized 
disk, generate Alfven waves and "feed" the disk corona 
with them (Varnier fc Tagger 2002). 

Currently sophisticated numerical codes are de- 
veloped, which allow to study the propagation of 
Alfven waves along an open magnetic flux tube 



(3aito et al. 2001). These simulations are aimed to clar- 
ify mechanisms of the coronal heating and of the for- 
mation of solar plasma flow patterns: viz. solar spicules, 
macrospicules and solar tornadoes. More general numeri- 
cal tools, like FINESSE and PHOENIX, aimed to study 
waves and instabiliti es in different kinds of flows, were very 
recently developed ( Bclien et al. 2002 ). We assume that 
all these numerical tools could also be used in a number 
of astrophysical applications involving "parent" flows and 
"inborn" Alfven waves interacting with each other. 

Most of above-cited studies considered simple (plane- 
parallel or locally plane-parallel) kinds of flows and the 
processes were treated by means of usual normal-mode 
analysis. The purpose of this paper (paper I) is to study 
'nonmodal' evolution of Alfven waves in swirling flows and 
to show that these flows may operate as efficient "ampli- 
fiers" of Alfven waves. 

In the next section we develop the mathematical for- 
malism and derive general equations governing the evolu- 
tion of Alfven waves in helical flows. However, the third 
section is dedicated to the study of the simpler example 
of a parallel SF ('pure outflow"). Still, even in this sim- 
ple case, we find that the amplification of Alfven waves 
takes place. The amplification mechanism is linear, tran- 
sient and can be described by appealingly simple mathe- 
matics. The amplification occurs within a relatively brief 
time interval and it appears as an abrupt, burst-like in- 
crease of the wave amplitude. It is tempting to argue that 
individual acts of wave amplifications, occurring sequen- 
tially, may ignite the "chain reaction" of nonmodal cascade 
amplification of Alfven waves. 

The fourth section of the paper is dedicated to fully 
helical SF. We show that when SF are helical ("outflow 
-t- rotation") there appear new classes of shear instabil- 
ities capable of generating high- amplitude Alfven waves. 
The presence of differential rotation is crucial for these 
instabilities and in certain cases the instability is of para- 
metric nature. These instabilities are powerful, often have 
a 'resonant' nature and are exclusively related with the 
kinematic complexity of 'parent' shear flows. 

The fifth section of the paper contains discussion of the 
obtained results and of their possible astrophysical appli- 
cations. We argue that such wave amplification processes 
may have various astrophysical manifestations including 
large amplitude Alfven waves actually observed in the so- 
lar atmosphere (Balogh et al. 1995). They could provide 



with these fiows (Ryutova et al. 2001). In accretion disks 
Alfven waves also seem to be actively interacting with 
the disk fiow - they are found to be unstable both in high- 



the necessary seeding for the development of MHD turbu- 
lence in hydromagnetic shear fiows. These processes may 
be present and may lead to perceptible morphological vari- 
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ety and diverse observational appearances in various kinds 
of 'accretion- ejection' flows: innermost regions of accre- 
tion disks, disk-jet transition regions, accretion columns 
in X-ray pulsars and cataclysmic variables, inner regions 
of galactic gaseous disks, etc. 

This paper will be followed by the second one (paper 
II), where we will consider the same flow structure but 
allow the perturbations to be fully compressible, bringing 
on-stage another two additional linear MHD wave modes 
- the slow and fast magnetosonic waves. 



2. Main consideration 

Our aim is to develop a theory of collective phenomena in 
helical flows of magnetized plasmas. We adopt the stan- 
dard MHD model and consider incompressible perturba- 
tions in an axisymmetric, cylindrical, steady flow of a 
plasma with uniform density (po = const), embedded in 
a vertical, homogeneous Bo = (0, 0, Bq) magnetic field. 
With p ^ po + Q, P ^ Po+P, V = U+u, B = Bo + B', the 
basic set of MHD equations for linearized perturbations is: 



Vu = 0, 



(1) 



Au+ (u- V)U = -— Vp- ^^x(VxB'), (2) 
Po 47r/9o 



AB = (B' • V)U + (Bo • V)u, 



VB 



0, 



(3) 
(4) 



where Dt=dt + (U • V) is the convective derivative oper- 
ator and U(r) is the vector field of the steady state flow. 
Unfortunately the actual kinematic portrait of astrophys- 
ical helical flows is largely unknown. There is no credible 
data for recently discovered solar macrospicules with heli- 
cal plasma motion {solar tornados) , for accretion columns 
related to X-ray pulsars and cataclysmic variables, for 
accretion- ejection flows and galactic and extragalactic jets. 
Therefore, we need to adopt a phenomenological model for 
U, that has to be general enough to encompass different 
possible sorts of real helical flows. Preferably the model 
must be three-dimensional and it should include a possibil- 
ity of both outflowing/inflowing and rotational modes o f 
motion. The model we choose is ( Rogava &: Poedts 1999|) : 



U(r)=[0, rn{r), U{r)], 



(5) 



where r = {x^ +2/^)^^^ is a distance from the rotation axis. 
For the angular velocity we take: 



n{r)^A/r'\ 



(6) 



with A and n as constants. This model implies the follow- 
ing two limiting cases: 

— Rigid rotation, with n — Q and ^ = f2o- 

— Keplerian rotation, with n = 3/2 and A = (GM)^/^. 



The corresponding Cartesian components of the linear 
azimuthal velocity = rfi(r) are: 



U,{x,y)^^Ay/r^ 
Uy{x,y) = Ax/r"". 



(7a) 
(76) 



A basic role in our analysis is played by the shear ma- 
trix Sik=dUi/dxk. Evidently rotational part of the flow 
velocity generates the following four nonzero components 
of the shear matrix: 



cr=Sxx = 2Anxoyo/rQ 



.n-i-2 



^yy — 



Al=Sxy - 

A2 = Syx 



-A[xl + {l-n)yl]/r-^+\ 
A[{l-n)xl + yl]hr\ 



(8a) 

(86) 
(8c) 

m 



As regards the outflow component of the velocity we 
can use the model ( Rogava et al. 2000 ): 



U,(x,y) = U^[l-{r/R)\ 



(9a) 



comprising a parabolically sheared "outflow" along the 
Z axis, being similar to the well-known Hagen-Poiseuille 
flow. The outflow model may be modeled also in a number 
of other wa ys. In astroph ysical jet literature, for instance, 
the model ( [Ferrari 1998| ): 



= U,r^/cosh[r" 



(96) 



is often used. One can easily see that the outflow compo- 
nent, given by (9), generates the following two components 
of the shear matrix: 



Ci=Sz 



^2UrnXQ/R^ 



C2=Szy — —2Umyo/R^, 



(10a) 
(106) 



Therefore, the resulting helical flow is the superposi- 
tion of a cylindrical outflow along the Z axis with a dif- 
ferential motion in the transverse cross-section of the jet. 
The complete 3x3 traceless shear matrix is: 



S = \ A2 -a 
Ci C2 , 



(11) 



This shear matrix has a number of interesting proper- 
ties. Its square is equal to: 



15^ 




(12) 



where we use the notation: 

£1=^2^2 + (tCi, (13a) 

e2=AiC7i -CTC2; (136) 
T={a^ +AiA2)^'^ = V^^n{r), (13c) 
and its cube is: H^^H = r^ll^ll. 
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The nonmodal local method for studying the dy- 
namics of linearized perturbations in kinematically com- 
plex flows (Lagnado et al. 1984, Craik & Criminale 1986, 
Mahajan & Rogava 199£) allows to reduce the initial set 
of partial differential equations for perturbation variables 
F{r,t), deflned in the real physical space, to the initial 
value problem formulated for the perturbation variable 
amplitudes F(k, f), defined in the space of wave numbers 
(k-space). The key element of this approach is the time 
variability of k's, imposed by the presence of the shear 
flow! This variability is governed by the following set of 
equations: 

k(i)+5^-k = 0, (14) 
henceforth for an arbitrary function / we use the notation: 

For the helical flow we find that k^^) = T^k^^). The 
corresponding characteristic equation dells'^ — AxX] = 
(with I being a unit matrix) yields the equation for the 
eigenvalues: 

A(A-r)(A + r) = 0, (15a) 



with the solutions: 



Al = 0, A2 



±r. 



(156) 



We see that the differential rotation parameter n plays 
a decisive role in determining the evolution scenario for the 
wave number vector k(t): when n < 1 (including the rigid 
rotation case) F is imaginary and the time evolution of 
k(t) is periodic, while when n > 1 (including the Keplerian 
rotation regime), F is real and makes the time behavior of 
k(t) exponential. 

In the dimensionless notation WK,a,z=Wx,j/,z, ^=*p/poi 
bx,y,z='iB^ y^^/ Bq, perturbation amplitudes evolve as: 



.(1) 



(k-v) = 0, (16) 

S-v^^kV -Cl[kb,^hh], (17) 

b(i) =5•b-fc3V,(^-:r,2;) (18) 

(k-b) = 0. (19) 

In the next two sections we separately consider the 
problem of a purely ejectional/injectional flow, and a flow 
with helical motion of plasma particles. We shall see that 
in both cases the Alfven waves sustained by the flow are 
amplifled. For the former case the amplification is tran- 
sient (algebraic instability), while for the latter case the 
velocity shear makes Alfven waves exponentially unstable. 

3. Ejectional flow 

This case corresponds to a — Ai = A2 = 0. We can in- 
troduce a vector C= [C^ = Ci,Cy = C2,Cz — 0], which is 
proportional to the position vector rg = (xo,j/O:0) con- 
necting the point A{xQ,yo, zq) (the centre of our local 
frame) with the Z axis in a z — zq plane. Its absolute 
value is equal to dUz/dr. The wave number dynamics is 
given by simple linear solutions: ki{t) = ki{0) — Cikzt, 



with i — X, y, while the vertical component kz stays con- 
stant. The symmetry of the flow implies that the vector 
product of C and k has a constant vertical component 
A= — (Cxk)2 = Cykx{t) — Cxky{t) = const(t), which, 
in its turn, is a linear combination of transverse compo- 
nents of the wave number vector. This fact suggests the 
decomposition fc||(i)=(C • k)/|C|, k±={C x k)2/|C|, i.e., 



k(i) = (fc||(t), kj 



(20) 



with the entire time dependence of k(t) contained in 

fc||(i) = fc||(0)- |C|M. (21) 

The system admits a surprisingly efficient and com- 
plete analytic description if we introduce, first, vectors 
of hydrodynamic (r2h=kxv) and magnetic (r2rn=kxb) 
vorticity and consider the following scalar products: 
V={C-v), V=(C-n,0, B={C-h), and B={C-nm)- The 
usefulness of these variables becomes apparent when we 
check that they obey the following set of first-order equa- 
tions: 



=-[?n(|k(i)|2)](i)T/ 
V^i) = AV^ + CikzB, 

= -kzV, 
^(1) = -AB-kzV, 



CikzB, 



(22a) 
(226) 
(23a) 
(236) 



which is remarkably simple with only one time-dependent 
coefficient (time dependence exclusively due to shear) in 
Eq. (22a). 

These variables give complete description of the sys- 
tem, because "physical" variables v and b are readily ex- 
pressed by means of V, V, B, and B via the following 
vector identities: 



1 



|Cxk|2 
1 

ICxkP 



v(Cxk) -Hy(kx(Cxk)) 



S(Cxk) + S(kx(Cxk)) 



(24a) 



(246) 



While for the kinetic and magnetic energies of the sys- 
tem we also have rather transparent expressions: 



1 



2|Cxk|2 



2|Cxk|2 

Etot=E]^ + 



B^ + B^M 



E„ 



(25a) 



(256) 
(25c) 



From (22-23) we can further notice that the Alfven 
waves, sustained by this system are "splitted" by the pres- 
ence of the velocity shear. 

— The first branch comprises just a simple, constant am- 
plitude and constant frequency Alfven mode, which is 
governed by the second order equation for B 



B^^'> + ojIB = 0, 



(26) 



with uJA=CAkz and with obvious solution B{t) ~ 
Csin{ujAt + (po)- 
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The second brancl£] is heavily modified by the pres- 
ence of the velocity shear. Its temporal evolution is 
described by the second order equation for B: 



Effetcive frequency square 



(27) 



We can further simplify this equation in terms of the 
variable 5'=|k|i3 which obeys 



(2) 



fc^lCxkl- 



= 0. 



Introducing the dimensionless parameter 



(28) 



(29) 



and the dimensionless "time" variable: 



(CA/|C|)fc||(t) = Ca[M - fc||(0)/|C|], (30) 



we may transform (28) into the remarkably simple dif- 
ferential equation: 



1 



= 0. 



(31) 



Solutions of this differential equation '^a{T) give the 
complete solution of the system, because all physical 
variables are readily expressed through this function, 
its first derivative and the trivial solution for the first, 
non-modified branch given by B{t). It means that the 
full solution of the initial value problem is expressible 
by certain combinations of these two Alfven modes: 
one that does not "feel" the velocity shear and the 
another one that is affected by it. 

Note that the symmetry of (31) implies that its so- 
lutions must be invariant with respect to the inversion 
operation. Namely, if a pair of initial values ^'a(^7o) 
and vE'i^^(— To) gives a certain solution then ^'^(To) and 
— ^i^-'(To) gives inversion-symmetric solution. 

Analysis of (31) begins by studying the temporal be- 
havior of the "effective frequency" 



f^cff(T)E 



(32) 



For sufficiently large absolute values of T, Qef / {T) ~ 
1, implying that at the corresponding periods of time, the 
system supports just the usual constant frequency and 
constant amplitude Alfven waves. However, when \a\ < 1 
there are two moments of time T±=±^ a(l — a) when the 
"effective frequency of shear-modified Alfven waves" be- 
comes zero and in the interval T_ < T < T+ it stays imag- 
inary. Therefore, within this time interval (with the width 
determined by Ar=2T+), swift and sudden changes in the 
temporal evolution of Alfven waves may be expected. 

^ This branch normally (i.e. in the absence of shear) also 
possesses constant Alfven frequency, but in a certain sense it 
is the incompressible limit of the fast magnetosonic wave. 



-2 -1.5 -1 -0.5 0.5 1 1.5 2 

T 

Fig. 1. The plots of the function il^fj (T) for different val- 
ues of the a parameter. Solid line corresponds to the case 
a — 0.9, dashed line to a = 0.7, dotted one to a = 0.5, 
and dashed-dotted to a = 0.3. 



Quantitative picture of this behavior is illustrated by 
Fig.l where we plotted the f^cff (^) different (but all 
\a\ < 1) values of a. One sees that with the decrease of 
a the minimum of ri^g(T) becomes sharper: the depth 
of the minimum in Fig.l (which accounts for maximum 
increments of the transient instability) steadily increases 
with the decreasing a. The width of the time interval, in 
which this function stays negative (and, correspondingly, 
n^g(T) is imaginary) is maximum for a = 1/2 and tends 
to zero when a~l and a ^ 1. 

At a first glance, from (31), it seems that decreasing a 
we would have the continuous increase of the transient 
amplification rate. However decreasing a we also make 
smaller the scaling factor between the physical time t 
and the variable T, which appears in (31). Besides, the 
small values of a imply smaller values of the Alfven speed 
and the vertical component of the wavenumber vector 
kz- Therefore, we can conclude that transient amplifica- 
tion factor for lower frequency Alfven waves is higher, 
but for the amplification to occur the system needs a 
longer time interval. This means, in turn, that outflows 
with shorter/longer lifetime values are expected to amplify 
higher/lower frequency Alfven waves. 

The presence of the imaginary effective frequency for a 
limited time interval means that Alfven waves excited and 
maintained by the shear flow become subject to a certain, 
velocity shear induced instability, which is present only 
within the limited time interval. That is why we specify 
this phenomenon by the term "transient instability. " Since 
the time interval is rather brief one can expect that the 
appearance of this instability will have a burst-like, explo- 
sive nature: initially at times T < T_ waves stay almost 
unaffected by the presence of the shear flow. However, 
as soon as the system will enter the transient instability 
interval T_ < T < T+, the Alfven waves undergo dras- 
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Psi, a = 0.1 




-100 -80 -60 -40 -20 20 40 60 80 100 



Fig. 2. The numerical solution of Eq.(31) featuring the 
function \E'o(T) for a = 0.1 case. The initial values are: 
*a(-100) ^ 0.01 and ^i^^lOO) = 0. 



Psi, a = 0.1 




vAAAAAAAAAAAAAA/ 



60 -40 -20 20 40 60 80 100 



Fig. 3. The numerical solution of Eq.(31), which is 
inversion-symmetric to the one given on Fig. 2. Initial val- 
ues here are: *a(-100) = -0.1354 and *1^^(-100) = 
-0.2186. 



tic and abrupt change in their amplitudes. Depending on 
the initial values, the mode of evolution will be a cer- 
tain mixture of swiftly decaying and/or increasing modes 
within the AT instability domain. But as soon as T > Tj^ 
the waves become stable again and their amplitudes do 
not change anymore The resulting wave amplitude is en- 
hanced/diminished in comparison with the initial ampli- 
tude depending whether transiently increasing/decreasing 
component was dominant for the initially excited wave. 

Numerical results, represented by Figs. 2-5, fully con- 
firm these qualitative expectations. Figs. 2 and 3 illus- 
trate the inversion symmetry of the functions 'I'a(T). The 
first of these two plots shows the temporal evolution for 
o = 0.1 and with initial conditions (—100) — 0.01 and 



^'i^^(-lOO) = 0. The figure is plotted for the interval 
— 100 < T < 100. Evidently the inversion symmetry im- 
plies that an another solution of this equation for the same 
value of a but with the initial values 'I'a(— 100) — \E'a(100) 
and 'i'i^'(-lOO) = -*i^\l00) will be exactly inverse- 
symmetric. In this particular example these initial val- 
ues are ^-al-lOO) = -0.1354 and ^'i^^(-lOO) = -0.2186. 
The inversion symmetry of these solutions is apparent. 
Physically this fact implies that the presence of the shear 
flow ensures burst-like and robust increase of amplitudes 
(energy) of some Alfven waves, while there are always 
other waves which, on the contrary, sharply loose their 
energy under the influence of the shear flow. 

Momentary appearance of the transient instability on 
the presented graphs is related to the narrowness of the 
transient instability interval, which is apparent from Fig.l. 
Note that the plotting time interval on Fig.l is taken very 
narrow in order to give magnified portrait of the behavior 
of ilcff, while Figs. 2-3 are deliberately drawn for the much 
wider range in order to illustrate the behavior of waves on 
a large time span. 

The family of solutions VPaCT) of (31) does not rep- 
resent a physical variable of the problem, so in order to 
recover information about the temporal evolution of per- 
turbations for physical variables it is more convenient to 
solve numerically (22-23) and to recover components of 
vectors v and b from (24). Besides, it is instructive to 
calculate the kinetic energy, the magnetic energy and the 
total energy of perturbations as given by (25). In order to 
track the temporal evolution of the shear-modified wave 
branch as such we set the amplitude of the unmodified 
component B(€) to be zero C = 0. Alternatively our nu- 
merical task was, first, to solve Eq.(28) in order to get 
functions \E'(<) and \['(^^(i). Second step was to calculate 
the physical variables by means of Eqs. (22-25). 

This set of calculations was performed for different val- 
ues of the system parameters and some representative ex- 
amples are given in Figs. 4-5. Note that they are plotted as 
functions of the real time variable t and not the variable T 
used in (31). From these figures we readily see that when 
A = (/cj_ = 0), the shear flow efficiently "pumps" energy 
into the longitudinal components of the velocity and the 
magnetic field perturbations, while the transverse com- 
ponents stay basically unchanged aside from the transi- 
tory, "burst-like" increase of their amplitudes in the brief, 
transient amplification phase. This is another indication 
of the above-mentioned fact that the velocity shear pri- 
marily affects the incompressible limit of the fast magne- 
tosonic waves. The asymptotic increase of the total energy, 
as shown in Fig. 4, being entirely due to the increase in Vz 
and fez, is quite substantial (about two orders of magni- 
tude for the given example). 

This is a typical mode of behavior for perturba- 
tions with A = and it is quite similar with the 
behavior of hydromagnetic waves in plane shear flows 
( Chagelishvili et al. 19931 ). By imparting a small but 
nonzero A (see Fig. 5), we can make all perturbation 
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Fig. 4. The numerical solution of Eq.(28) presented for 
all components of v (solid lines on first three plots) and b 
(dashed lines on first three plots) vectors. Last plot shows 
time evolution of Etot{t)/Etot{0)- The values of parame- 
ters are: Ca = 0.2, = 20, ky = 10, = 0.1, = 0.1, 
Cy = 0.05, a = 0.1789. Initial values are: *(0) = 0.01 and 
*(i)(0) = 0. 



Vy and by 






50 - 



Fig. 5. The numerical solution of Eq. (22-23) with the 
same parameters as on Fig. 4 except ky = 10.1. 

components grow. However, the overall increase of the to- 
tal energy of perturbations in the latter case is somewhat 
smaller than in the former case. 



4. Helical flow 

For more complicated MHD flows, the time evolution of 
the wavenumber vector k(t), governed by equations (14), 
becomes nonlinear and it makes the temporal behavior 
of linear perturbations much more complex. A relatively 
simple, 2-D hydrodynamic system of such complexity was 



recently investigated (Mahajan & Rogava 1999). For the 



considered in this section, we may deduce from (15) and 
(13c) that the perturbations may grow either exponen- 
tially (when differential rotation rate n > 1), or vary in 
time in a periodic way (when n < 1). Below we develop 
the mathematical formalism for the study of perturbation 
evolution for these two classes of solutions. 

4.1. General formalism 

The starting set of equations is again (16-19). Taking time 
derivative of (16) we find out that 



(k-v) 



(1) - _ 



(5^.k).v - k.(5-v) - |k|2(P + C'ih) = 0, 



which after noticing that (iS"^-k)-v = k-(5-v) implies that 
the following algebraic relation holds: 

P + Clb, = -2k-(5-V)/|k|2 = 2k(i)-v/|kp. (33) 

Note that this relation follows from (16), which, in 
turn, is direct consequence of the 'incompressibility' con- 
dition. The relation (33), applied to (17), helps to reduce 
the initial system to the following pair of the first-order 
ODE's: 

v(i) + 5-v = Cifc.b - 2(k(i)-v)k/|kp, (34) 

b(i) =5-b-fc^v. (35) 

It is worthwhile also to note that k(k'^^'-v) = 
vx(kxk'^-'), which allows to rewrite (34) in the follow- 
ing form: 



.(1) 



Sv = C^/c^b + Rxv, 



with R defined as: 



R(t)=(kxk(i))/|k|- 



(36) 



(37) 



The interesting property of this form of the Eq. (34) is that 
the only time-dependent coefhcient is just represented by 
one vector R{t)l 

We can reduce (34-35) to the following explicit second- 
order equation for the magnetic field: 

b(2) + \col - 5^1 b + fk(i) •b(i) + k(2) -bl = 0. (38) 
|kp L J 

This equation is quite informative, because it shows 
explicitly the effects of the helical flow imposed upon the 
evolution of perturbations. Note that the term containing 
is also present in the pure outflow case (see for 
comparison Eq. (27)). While terms containing S^h and 
k(2) . b are characteristic to the case of helical flow and 
they reflect different aspects of shear-induced processes in 
helical flows. 

Obviously (38) is a general equation and we can base 
our analysis on it. However a different approach may also 
be useful. Let us introduce "projection variables" : 



MHD helical flow, specified by the shear matrix (11) and 



U —Sij ki Vj 



H=Sijkibj 



-k(i).v, 
-k(i).b, 



(39a) 
(396) 
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(39c) 
(39d) 



Using these definitions and our basic set of equations 
we find: 



|k|2 



B^^^ = -KV. 



(40) 

(41) 
(42) 
(43) 



From this set we can derive the following two second- 
order equations for H and B: 



H 



(2) 



ikp 



[H = -2B^'\ 



where 



Moreover, if we define 

we can rewrite (44) also as: 

(k-k(2)) 



(2) 



m 



^2|k|B 



(1) 



(44) 
(45) 

(46) 
(47) 

(48) 



Note the appearance of the vector R in this equation. 
Overall, (45) and (48) give rather transparent analytic pre- 
sentation, which can be used for the analysis of interesting 
classes of solutions. Note also that in the absence of rota- 
tion f{t) = these equations become decoupled and they 
reduce to the "pure outfiow" case. In that case, we had 
a splitting of the mode into the usual Alfven and shear- 
modified ("hybrid" as we can also call it) Alfven mode. 
Below we shall see that rotation not only gives birth to 
new instabilities but it couples these two branches of the 
Alfven mode! 

4.2. Exponential k(t) case 

Depending on the actual type of the motion in the plane 
normal to Z there will be different regimes of the time 
variation of k(i) (see for details Mahajan fc Rogava 19991 ). 



Instabilities are present both when the wavenumber varies 
exponentially and periodically. 

In the former case, i.e., when the absolute value of 
k(i) increases exponentially, the dissipation effects will 
eventually start to be important and try to damp the 



mode (Mahajan & Rogava 1999) despite the presence of 
the shear instability. The same argument is valid for those 
cases, too, when the temporal growth of the wavenumber 
is linear. Therefore, considering shear-induced nonmodal 
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Fig. 6. The numerical solution for the case when temporal 
evolution of k's is exponential. The values of parameters 



are: kx ■ 
A2 = 10- 



5, hy 



10, 



1, (J 



Ci = 0.8, C2 = 0.9, Ca 



= 10-2, Al 
0.5. 



0.3, 



effects in inviscid fluids in cases when |k| is monotonously 
increasing, we should realize that the description is phys- 
ically meaningful only for initial (finite) times. 

It is worthwhile to make a look at some of these solu- 
tions. On Fig. 6 we present the results of numerical calcu- 
lations, taken for the case, when kx — 5, ky — 10, kz — I, 
(7 = 10-2, Al = 0.3, A2 = 10-2, Ci = 0.8, C2 = 0.9, 
Ca — 0.5. The value of F^ = 3.1x10-^ is positive and the 
value of |k(i)| is exponentially increasing (see Fig. 6a). The 
characteristic amplification time scale T~Vt~^ is of the or- 
der of the period of rotation. The resulting wave evolves 
accordingly, which is clearly visible from the figures. It 
is noteworthy to see that, analogous to the pure outflow 
case, shear flow energy is most efficiently absorbed by the 
longitudinal components of the velocity and the magnetic 
field perturbations. Still, now, the energy growth, plotted 
again for the -Etot(i)/-E^tot(0), is not algebraic but expo- 
nential and is entirely due to the exponential evolution of 
the perturbation wave number vector. 

4.3. Periodic k(t) case 

When values of F^ are negative, the temporal evolution 
of k(t)'s is periodic, which means that the absolute value 
of the wavenumber vector stays bounded within certain 
limits. It means that in this case shear-induced effects 
have potentially much wider durability and may poten- 
tially lead to strongly perceptible effects. The limits of 
this paper does not allow us to present full analysis of 
all possible classes of solutions. The most interesting fea- 
ture of this case is that Alfven waves become paramet- 
rically unstable in a helical flow. The Figs. 7-9 demon- 
strate this noteworthy fact. All figures are drawn for the 
case when Ai = —A2 — 0.3, and a = 0.01, so that 
VF=(-F2)1/2^0.2998. Fig.7 shows the case when Alfven 
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Fig. 7. The numerical solution for the case when Ai = Fig. 9. The numerical solution for the case when 
-A2 = 0.3, a = 0.01, Ci = 0.8, C2 = 0.9, W~0.2998 and VF~0.2998 (the same set of parameters as on the Fig.7) 
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Fig. 8. The numerical solution for the case when 
W^~0.2998 (the same set of parameters as on the Fig.7) 
but with u)A = 0.3. 



frequency is oja = 0.29. The system is still out of the para- 
metric resonance and the system, therefore, exhibits just 
Alfven waves, modulated by the periodic shear. The pic- 
ture is drastically different when wa = 0.3 — the Fig. 8 
shows that there appears a robust instability of paramet- 
ric (resonant) nature that leads to the strong exponential 
increase of the system energy in the same span of time as 
on Fig.7. The resonance is quite sharp, because as soon as 
Alfven frequency slightly exceeds the value of W the sys- 
tem becomes stable again, as it is visible from the Fig. 9, 
which is drawn for the case when lua = 0.31. 

Further study shows that actually there are several 
isolated regions of instability, which happen to be centered 
around the values coa = nW, n = 1,2,.. and usually first 
resonance W = wa is the most efficient one. This is a 
strong indication in favor of the resonant and parametric 



Fig. 10. The numerical solution for total energy of per- 
turbations normalized on their initial values. The values 
of parameters are: = 20, ky = 5, kz = 1, a = 0.6, 
Al = 0.5, = -0.9, = 0.01, Cy = 0.05. In this case 
W = 0.3, while wa = 0.3,0.5,0.6, and 0.8 on the plots 
labeled (a), (b), (c), and (d), respectively. 



nature of this instability. Fig. 10 and Fig. 11 illustrate this 
interesting property of the system. The plots are drawn for 
the total energy of a perturbation, normalized on its initial 
value: ii^tot(i)/-^tot(0). The simulation is made for the case 
when = 20, ky = 5, fc^ = 1, cr = 0.6, Ai = 0.5, A2 = 
-0.9, Ca, = 0.01, Cy = 0.05. Note that in this case W = 
0.3. Making simulations for different values of the Alfven 
velocity we see that perturbations are strongly unstable 
when LUA—nW. Namely the instability is present on Fig. 10 
(a) and (c) and Fig. 11 (a) and (c), where the values of wa 
are 0.3, 0.6, 0.9, and 1.2. respectively. At the same time we 
see that for intermediate values of ua, when nW < wa < 
{n + l)W, the Alfven waves stay stable. Unfortunately 
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Fig. 11. The numerical solution for the total energy of 
perturbations. The values of parameters are the same as 
on Fig. 10. while — 0.9,1.1,1.2, and 1.4 on the plots 
labeled (a), (b), (c), and (d), respectively. 

the complexity of the system does not allow us to perform 
strict analytic analysis and locate the actual width and 
size of instability regions. However, our results do show 
that there are several such regions and presumably the 
efRciency of the instability decreases with the increase of 
n. 

5. Discussion 

The results of the present study clearly show that both 
helical and purely cjectional flows of MHD plasmas could 
serve as efficient "amplifiers" of Alfven waves. In alliance 
with other shear-induced phenomena (such as shear- 
induced wave transformations and/or generation of shear 
vortices), these processes could lead to the cascade ampli- 
fication of hydroma gnetic waves leadin g to MHD turbu- 
lence in these flows ( Rogava et al. 2000|) . In more general 
terms it could be stated that the presence of the velocity 
shear strongly modifies the nature of Alfven waves such 
that they can extract energy from the mean flow. 

Certainly, while talking about these processes, we 
should bear in mind that due to the limitation of the 
'nonmodal approach' we are not able to describe quan- 
titatively the onset of these processes in the real physical 
space. Besides the linearized nature of our calculations 
prevents us from taking into account of the back reaction 
of shear-modified Alfven waves on the mean flow. The last 
but not the least limitation is related with the neglect of 
the dissipation in the flow. 

Can Alfvcnic perturbations still extract energy from 
the background flow when dissipation is present? For 
magnetoacoustic modes sustained by a simple, linearly 
sheared two-dimensional flow it was recently found 



(Bodo et al. 2001) that the presence of dissipation doesn't 
prevent a shear-modified MHD mode from the energy 
extraction. For the case of helical dissipative flows until 



now no numerical studies were done. Numerical study of 
these processes in the swirling flows requires 3-dimensional 
MHD simulations over large time intervals, which is quite 
a tough challenge. We expect that the outcome of the 
"competition" between the increase of the energy due to 
the presence of the shear flow and viscous damping due 
to the presence of the dissipation will depend on the time 
scales of both processes and initial conditions. Relevant 
numerical studies are planned to be started in the near 
future. 

However, basing on our results about 'nonmodal' tem- 
poral evolution of Alfvenic perturbations, we can give 
qualitative picture of the physical appearance of these 
Alfven waves in the presence of viscous dissipation for all 
those three subclasses of flows - viz: (a) ejectional outflows 
with no rotation; (b) helical flows with exponentially in- 
creasing |k(t)|; (c) helical flows with periodically changing 
|k(i)|) - which were considered. 

For an ejectional outflow, with linear temporal growth 
of the |k(i)|, a monochromatic Alfvenic disturbance would 
undergo one act of transient, 'burst-like' amplification fol- 
lowed by the viscous decay. For the package, containing 
harmonics of different |k(0)|, one could expect a sequence 
('cascade') of their transient amplifications mixed with 
their viscous damping. This sort of the Alfven wave, while 
passing through the flow, would slowly decay but repeat- 
edly regain its energy via the chain-sequence of amplifica- 
tion events. 

For helical flows, where the temporal increase of the 
|k(t)| is exponential, the Alfven waves could be less for- 
tunate. It is true that their energy might grow exponen- 
tially, but the increase of the |k(t)| seems to be a more 
rapid process than the increase of the total energy of the 
Alfven wave (see Fig. 6). Therefore, even if initially the 
Alfven waves could extract some energy from the back- 
ground flow, the powerful viscous damping would sooner 
or later prevail and return the energy back to the flow in 
the form of heat. It means that in this situation, perhaps, 
one could not see any high amplitude Alfven waves, but 
one would rather witness an efficient heating of the flow 
due to this 'self -heating' process: when Alfven waves serve 
as energy transmitters, extract a part of the flow's regular 
kinetic energy and returning the energy back to the flow 
in the form of the heat. 

For helical flows, where the variability of the |k(t)| 
is periodic, the situation can be the most exotic. The 
strength of the viscous damping in such flows is limited, 
because the value of the |k(t)| is kept finite: it doesn't in- 
crease neither algebraically nor exponentially. At the same 
time, Alfven waves are modulated by the presence of the 
flow and there exist separate regions of parametric insta- 
bility, which implies that under favourable conditions the 
Alfven waves propagating through this kind of flow can 
be unstable. Therefore, flow patterns of this gender might 
serve as efficient amplifiers of Alfven waves. 

In all these cases the ultimate 'fate' of any individual 
Alfven wave would depend on the time- and length-scales 
of their propagation and the flow geometry, as well as, of 
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course, on the actual strength of the viscous dissipation in 
the flow. 

The shear-modified Alfven waves seem to be the hy- 
brids of the usual Alfven waves and the so called 'Kelvin 



modes' ( [Marcus fc Press 1977| , [Rogava et al. 1998[ ). The 
latter are nonperiodic, transiently growing modes inherent 
to the incompressible neutral shear flows. In the absence 
of the equilibrium magnetic fleld, the Kelvin Transients 
are the only modes of collective behaviour exhibited by 
the flow. The shear-modified Alfven waves are formed by 
combining the features of the usual waves with the tran- 
sient shear vortices; their ability to extract/give energy 
from/to the mean flow follows from the latter. 

Within the framework of the nonmodal approach we 
are able to track and study properties of these shear- 
modifled Alfven waves in the space of wave numbers (k- 
space). The orientation and magnitude of a wave num- 
ber vector, being variable in time and governed by Eq. 
(14), varies because of the local straining of the back- 
ground flow. The evolution of the k(t) fleld, in its turn, 
exercises influence on the temporal behaviour of physical 
variables: density and pressure perturbations, components 
of the magnetic fleld and the velocity vectors. In terms 
of the usual physical coordinates (space and time) these 
variables are periodic in each space coordinate. However, 
these modes are 'non- normal' ( Butler fc FarrcU 1992[ ) be- 
cause both wavenumbers and amplitudes associated with 
each disturbance are functions of time, implying that the 
solutions are not separable in space and time. The lat- 
ter, nonexponential time-dependence of the solutions is 
disclosed through the nonmodal approach and it helps 
to "smoke-out" new classes of modes hardly accessible 
through traditional normal mode analysis. In order to 
have complete real-space physical picture of the evolu- 
tion of these Alfvenic perturbations we need, therefore, 
to develop phcnomcnological approach and try to un- 
derstand the dynamics in terms of vorticity dynamics 
( Craik fc Criminale 1986 ). Alternatively one might ap- 
proach the problem through direct numerical simulations, 
generalizing the very first attempts of numerical studies of 
much simpler MHD velocity patterns ( Bodo et al. 2001 ). 

Turning our attention back to astrophysics, we might 
argue that the role of the shear-modified hybrid Alfven 
waves could be quite important in a number of astro- 
physical situations. Let us first assess the relevance these 
waves could have to the solar plasma flows. It is known 
that higher layers of the solar chromosphere consist of 
long, vertical col umns (spicu les) rising above the gen- 
eral background (Athay 1986) and well- visible at the ex- 
treme limb of the Sun, where the overall structure of the 
chromosphere resembles a "burning prairie." The spicules 
rise vertically, like slender magnetic fiux tubes (filaments) 
out of the network beginning in the chromosphere and 
threading through the solar transition region into the low 
corona. Their average lifetime is 10-15 minutes, character- 
istic length scale lO.OOOfcm, and characteristic diameter 
700— lOOO/cTO. Often the spicules appear to have an "Eifel 
tower" shape with vertical flow velocities in the range 



20 — 25km/ s. There is also some evidence for plasma ro- 
tation in spicules (Kudoh fc Shibata 1999). It means that 
the plasma motion within spicules may easily be of heli- 
cal nature. The recent discovery of swirling macrospicules, 
called solar tornados (Pike & Mason 1998) provides addi- 
tional evidence in this direction. 

If we imagine flow patterns, like spicules or solar tor- 
nadoes, with moderate storage of total (quasi)equilibrium 
energy it may happen that transiently amplifled Alfven 
waves, bred within these structures, could efflciently 'suck 
out' much of their energy and lead to eventual exhaus- 
tion and disappearance of the flow patterns. The shear- 
induced cascade of Alfven waves could take away much 
of the available energy of the "parent" flow. The observa- 
tional consequences may be dramatic; we could imagine 
individual spicules acting as short-lived plasma cannons 
firing away heavy blasts of large-amplitude Alfven waves 
and disappearing after the wave package has drained most 
of the energy from the "parent" flow. 

It is quite fascinating to surmise that cascade of tran- 
sient ampliflcations of Alfven waves within spicules, to- 
gether with other kinds of shear-induced wave processes 
(like, e.g., "MHD wave oscillations" Rogava et al. 2000) 
may account for the quasisteady appearance of spicules 
in the solar atmosphere. Taking the "burning prairie" 
metaphor a little further, we may speculate that shear 
induced processes, "burning" individual "grass blades" , 
could be responsible for the overall conflagration of the 
spicule network which, in turn, provides an efficient mech- 
anism for the transfer of energy from the quasisteady mo- 
tions of solar plasmas into the fiux of large-amplitude 
Alfven waves. These waves are able to penetrate through 
the solar transition region and reach the lower layers of 
the solar corona. Here the processes of the she ar-induced 
"MHD wave oscillations" ( Rogava et al. 200C ) may con- 
tribute to the appearance of compressible MHD modes 
(slow and/or fast magnetosonic waves) which, in turn, are 
effectively damped giving back their energy to the solar 
plasma particles, accelerating them and providing initial 
launch of the solar wind flows. This highly speculative 
picture must be carefully developed and tested in order 
to make it convincing. We are in the process of develop- 
ing a numerical model for an individual helical solar flow 
pattern, and hope to see the above-described processes in 
real physical space through numerical simulations. 

The parametric Alfven instability may be effective in 
flows, where conditions are favorable for the correlation 
between the rate of rotation and stretching of the flow 
lines (which determines the value of the characteristic fre- 
quency W) and the strength of the equilibrium magnetic 
field (which determines the magnitude of the Alfven fre- 
quency llJa)- Such a correlation may easily happen within 
solar jet fiows, where both the ejectional and rotational 
modes of motion are present. Any time such correlation 
takes place it will lead to a burst-like generation of high- 
amplitude Alfven waves with corresponding "exhaustion" 
(and, maybe, even disappearance) of the "parent" helical 
flow pattern. Since the modern apparatus used on the last 
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generation solar satellites enables one to measure both the 
strength of magnetic fields and the rate of rotation, one 
can assume that it will be fairly realistic to verify whether 
the appearance of high-amplitude Alfven waves is really 
related to shear-induced resonant parametric instability of 
helical flows. 

We think that shear-induced instabilities of this na- 
ture may be also present in accretion- ejection flows. The 
differential rotation parameter n controls the k(i) dynam- 
ics which, in turn, specifies the nature of fluctuations 
in a given flow. When n < 1 (including the rigid rota- 
tion case) r is imaginary and k(t) is periodic. In such 
cases (e.g., innermost regions of galactic gaseous disk, be- 
lieved to have almost constant rotation rates) the system 
may sustain parametrically unstable Alfven waves. While 
when n > 1 (including the Keplerian rotation regime), 
r is real and makes the time behavior of k(i) exponen- 
tial. This regime can be realized in different kinds of 
quasi Keplerian accretion disks. Note that the charac- 
teristic amplification time scale T~i7~^ is of the order 
of the inverse period of rotation, being at least of the 
same order as the "Magnetor otational instability" growth 
rate ( Balbus fc Hawley 1998 ). Since the latter instability 
is thought to be the strongest one, accounting for the 
turbulence in accretion disks, we could surmise that the 
shear-induced instability discussed in this paper could also 
play a role in the onset of turbulence in accretion-ejection 
flows. 

Finally, considering our results in the context of galac- 
tic gaseous disks and their large-scale magnetic fields, 
we might envisage a possible relation between the shear- 
induced amplification of Alfven waves and standard dy- 
namo action in swirling flows. Within galactic disks al- 
gebraic and/or usual instabilities of the nonmodal origin 
might provide the important mechanism for the initial am- 
plification of the magnetic energy and this process could 
compete with the omnipresent Ohmic diffusion. This is 
a challenging task to see what is the role of nonmodal, 
shear-induced processes in the framework of the general 
theory of dynamo generation of galactic magnetic fields. 
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